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Abstract 

The strength functions of isovector dipole resonances in the even oxygen 
isotopes 18 ~ 24 are calculated by using an extended version of the time- 
dependent Hartree-Fock theory known as the time-dependent density-matrix 
theory (TDDM). The results are compared with recent experimental data 
and also with a shell-model calculation. It is found that the observed isotope 
dependence of low-lying dipole strength is reproduced in TDDM when the 
strength of the residual interaction is appropriately chosen. It is also found 
that there is a difference between the TDDM prediction and the shell-model 
calculation for 24 O. 
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The study of structures and reactions of nuclei far from stability has been one of the 
most active fields of nuclear physics in the past decade jj]]. The microscopic description of 
such nuclei will lead to a better understanding of the interplay among the strong, Coulomb, 
and the weak interactions as well as the enhanced correlations present in these many-body 
systems. For neutron rich nuclei, in the extreme limit, the proton and neutron densities may 
exhibit strikingly different behavior. One of the observables for studying this phenomenon is 
the evolution of the isovector dipole modes along an isotope chain. Several years ago explo- 
rative calculations were performed for isovector dipole (El) and isoscalar quadrupole (E2) 
resonances in the unstable oxygen isotope 22 based on the time-dependent density-matrix 
theory (TDDM) 0. It was found that there are low-energy El and E2 modes associated 
with excess neutrons. At that time there were no experimental data for comparison and 
the TDDM calculations were not highly quantitative since spin-orbit force was neglected in 
the calculation of a mean-field potential, and the neutron 2s and Id states were assumed to 
be fractionally occupied, apparently overestimating the effect of ground-state correlations. 
Recently the isotope dependence of low-lying El strength has been measured at GSI || 
for oxygen isotopes and compared with a shell-model calculation ||. The aim of this pa- 
per is to present more quantitative calculations of El modes for these isotopes using the 
time-dependent Hartree-Fock (TDHF) and the TDDM theories. Two recent improvements 
make it possible to perform better TDDM calculations: one is the treatment of ground-state 
correlations. As will be described below, an adiabatic treatment of the residual interaction 
enables us to obtain a correlated ground state which is a stationary solution of TDDM || . 
The other improvement is the extension of the TDDM program to include spin-orbit force, 
as it was done for TDHF calculations ||. 

TDDM is an extended version of TDHF and is formulated in order to determine the 
time evolution of one-body and two-body density matrices p and p<i in a self-consistent 
manner J7|. TDDM, therefore, includes the effects of both a mean-field potential and two- 
body correlations. The equations of motion for p and p 2 can be derived by truncating the 
well-known BBGKY hierarchy for reduced density matrices ||. To solve the equations of 
motion for p and pi , we expand p and C 2 (the correlated part of p 2 ) using a finite number 
of single-particle states ip a which satisfy a TDHF-like equation, 

p(ll',t) =J2n aa ,(t)-ip a (l,t)r a ,(l',t), (1) 

aa' 

C 2 (12l'2',t)=p 2 -A(pp) 

= J! Ca/3a'l3' (t) 
a/3a'(3' 

x Mht)M^t)r a ^',t)^(2',t), (2) 

where the numbers denote space, spin and isospin coordinates. Thus, the equations of 
motion of TDDM consist of the following three coupled equations : 

ih^ a (i,t) = h(i,t)*jj a (ht), (3) 

iM aa > = Y^[( a P\ v \l8)C lSo ,r/3 - C al3 ^5(lS\v\a'P)], (4) 

/3 7 <5 
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ihC a j3a'l3' — B a! 3 a if3> + P a j3a'(3> + Haj3a'f3', 



(5) 



where h is the mean-field Hamiltonian and v the residual interaction. The terms on the 
right-hand side of Eq.(5) contain all the two-body correlations including those induced by 
the Pauli exclusion principle 0. The TDDM equations of motion satisfy conservation laws 
of the total number of particles and the total momentum and energy. The small amplitude 
limit of TDDM was investigated |J and it was found that if only the 1 particle - 1 hole and 
1 hole - 1 particle elements of n aa > and the 2 particle - 2 hole and 2 hole - 2 particle elements 
of C a p a ipi are taken, the small amplitude limit of TDDM is equivalent to the conventional 
second RPA (SRPA) jlQj . Thus, TDDM is a more general framework than the conventional 
SRPA. 

The El strength function is calculated according to the following three steps: 
1) A static Hartree-Fock (HF) calculation is performed to obtain the initial ground state. 
The Skyrme III with spin-orbit force is used as the effective interaction. It has been reported 
n| that the Skyrme III is useful even for very neutron rich nuclei. Unoccupied single-particle 



states up to the 2s and Id orbits are also calculated for both protons and neutrons to solve 
the TDDM equations for n aa i and C a p a ipi. The single-particle wavefunctions are confined 
to a cylinder with length 20fm and radius lOfm. (Axial symmetry is imposed to calculate 
the single-particle wavefunctions ||.) The mesh size used is 0.5fm. The neutron 1^5/2 state 
is assumed to be partially occupied to obtain the HF ground states for 18 and 20 O. 

2) To obtain a correlated ground state, we evolve the HF ground state using the TDDM 
equations and the following time-dependent residual interaction of the ^-function form 

v (t) = v {l - e- t/T )5 3 (r - ?). (6) 

The time constant r should be sufficiently large to obtain a nearly stationary solution of the 
TDDM equations ||. We choose r to be 150fm/c. The strength of the residual interaction 
is determined to approximately reproduce the observed occupation probability of the proton 
1g?5/2 state in 16 The obtained value of Vq is — 230MeVfm 3 , which might be smaller than 
the value of approximately — 300MeVfm 3 found in the literature for pairing-gap calculations 
|T3|| . The time step size used to solve the TDDM equations is 0.75fm/c. The calculated 



occupation probabilities of the neutron 2si/2 and ld 3 / 2 are 2 ~ 3% in the oxygen isotopes 
considered here. 

3) The El mode is excited by boosting the single-particle wavefunctions at t = 5r with 
the dipole velocity field: 



M^r) — > e ikD ^ a , (7) 



where 



N Z 

D(z) = —z, ( — -z) for protons (for neutrons). (8) 

Here, Z and iV are the numbers of protons and neutrons, respectively, and A = Z + N. 
When the boosting parameter k is sufficiently small, the strength function defined by 

S(E) = J2\(^n\D\^o)\ 2 6(E~E n ) (9) 
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is obtained from the Fourier transformation of the time- dependent dipole moment D(t) as 

1 r°° Et 

s < B » = rts/o do) 

where 

D(t) = J D(z)p(f,t)d 3 f. (11) 

In Eq.(9) |$o) is the total ground-state wavefunction and |<&„) the wavefunction for an 
excited state with excitation energy E n . It is very time consuming to solve the TDDM 
equations (especially Eq.(5)) for a long period. Thus, we stop TDDM calculations at t — 
1200fm/c. The upper limit of the time integration in Eq.(10) is limited to 450fm/c. To 
reduce fluctuations in S(E), the dipole moment is multiplied by a damping factor e~ rt ^ 2h 
with T = lMeV before the time integration. Since the integration time is limited, the 
strength function in a very low energy region (E < 2irh/450 3MeV) is not well determined. 
The energy-weighted sum rule (EWSR) is expressed as 

J S(E)EdE = ^(%\[D, [H,D]]\$ ) 

h 2 ZN 
~ 2m A 

+ / p P (7)p n (7)d 3 r+ ]T J r a (m(^(f1M^d 3 fC a/ ^), (12) 

where m is the nucleon mass, and t\ and t 2 are the parameters for the momentum dependent 
parts of the Skyrme interaction. We assume that the Hamiltonian H consists of a two-body 
interaction of the Skyrme type. The first term on the right-hand side of the above equation 
corresponds to the classical Thomas-Reiche-Kuhn (TRK) sum rule and the second term, 
the enhancement term, is due to the momentum dependence of the Hamiltonian. The 
contribution of the momentum dependent part is about 28% of the total EWSR value in the 
oxygen isotopes considered here and the term proportional to C a /3 a >j3> describing the effects 
of ground-state correlations is quite small (less than 1% of the total sum rule value). 

In Figs.l~4, the El strength functions of 18 ~ 24 calculated in TDDM (solid line) are 
compared with those in TDHF (dotted line). The TDHF calculations presented here are 
equivalent to the RPA calculations without any truncation of unoccupied single-particles 
states because the TDHF equation for the boosted single-particle wavefunctions ip a is solved 
in coordinate space. The boundary condition for the continuum states, however, is not 
properly taken into account in our calculation because all the single-particle wave functions 
are confined to the cylinder. Therefore, the calculated strength functions slightly depend on 
the cylinder size. The difference between the TDDM and TDHF calculations is due to the 
effect of two-body correlations, which may have two aspects: one is to induce ground-state 
correlations which increase the occupation probabilities of weakly-bound neutron orbits. 
The increase in low-lying strength (E < 15MeV) seen in the TDDM results for all the 
isotopes is due to partial occupation of the neutron 2s 1 / 2 and ld 3 / 2 states. The comparison 
between the TDDM result for 24 and those for other isotopes suggests that the occupation 
of the neutron 2s\j 2 state is responsible for the increase in the El strength in the very 
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low energy region (around E = 5MeV). This is because the appearance of a prominent 
peak at 6MeV in 24 is related to nearly full occupation of the neutron 2s\/2 state. The 
other aspect of two-body correlations is to increase the width of the giant dipole resonance 
(GDR) which is located around 25MeV. There are background 2 particle - 2 hole states 
which consist of lfwo excitations of protons and Otuv excitations of neutrons. The coupling 
of the GDR to these states leads to the spreading of the GDR strength. The fractions of 
the EWSR values depleted over the energy range between OMeV to 40MeV are about 85% 
in TDDM and about 90% in TDHF. The EWSR values depleted below 15MeV are shown 
in Fig. 5. The ordinate indicates the ratio of the EWSR value to the TRK value. The 
results in TDDM with vq = — 230MeVfm 3 (rhombus) are smaller than the experimental 
data except for 22 0. In TDDM, the EWSR values in the low -energy region crucially depend 
on the strength of the residual interaction. Since the strength of the residual interaction is 
not well-known in unstable nuclei, we performed TDDM calculations using slightly stronger 
residual interaction with vq = — 330MeVfm 3 to increase El strength in the low-energy region. 
The value of v = — 330MeVfm 3 was used in our previous studies 
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reasonable spreading widths for giant resonances in stable nuclei |I4j]. As an example of 
the TDDM calculations with vq = — 330MeVfm 3 , we show the result for 20 O in Fig.6. The 
increase in the strength below lOMeV is due to the increase in the occupation probability 
of the neutron 2syz state and the enhancement of the peak at 15MeV is caused by further 
fragmentation of the GDR strength. The EWSR values depleted below 15MeV are shown 
by circles in Fig. 5. The magnitude of the low-lying strength is increased and now becomes 
close to the experimental data except for 22 0. The fact that the low- lying El strength in 
20 O is most increased with the use of the stronger residual interaction may be explained 
by the large fragmentation of the GDR strength as seen in Fig.6. The TDDM results for 
24 shown in Fig.5 significantly differ from the shell-model calculation 0. The difference 
originates in the fact that the shell-model calculation gives much smaller spreading width to 
the GDR in 24 than the TDDM and TDHF calculations do. Since the large amount of the 
El strength is concentrated in the GDR, the low-lying El strength calculated in the shell 
model becomes small in 24 0. 

Recently, surveying observed neutron separation energies and interaction cross sections, 
A. Ozawa et al. has pointed out that neutron number N = 16 becomes a new magic 
number for nuclei near the neutron drip line. Shell model calculations done by Y. Utsuno 
et al. [0 have also shown that a large energy gap between the 2s\/2 and ld^/2 orbits makes 
N = 16 a new magic number for neutron-drip line nuclei. They predict that the excitation 
energy of the first 2 + state in 24 becomes the largest among the unstable oxygen isotopes as 
a consequence of the large energy gap [TO. Since the lowest-lying El mode in 24 probably 
originates in the transition of a neutron from the 2sx/2 to the 2p^/2, its position and strength 
may also be used to explore the shell structure of 24 0. However, a quantitative discussion 
on the magicity of N = 16 is difficult for the present TDDM calculations. This is because 
the energy gap between the neutron 2si/ 2 and ld 3 /2 orbits does not significantly depend on 
N in our calculations. In order to explore the magicity of iV = 16 and its relation to the 
low-lying El modes, we might need larger single-particle space and more realistic residual 
interactions. 

In summary, the strength functions of the isovector dipole resonances in the even oxygen 
isotopes 18 ~ 24 were calculated by using an extended version of TDHF known as TDDM. 
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These calculations are much more improved and quantitative than the previous ones 
in two aspects: the inclusion of spin-orbit force in a mean-field potential and the use of a 
correlated ground state. The results were compared with recent experimental data and also 
with a shell-model calculation. It was found that if the strength of the residual interaction is 
appropriately chosen, TDDM approximately reproduces the observed low-lying El strength. 
It was also found that the TDDM prediction differs from the shell-model calculation for 24 0. 
It was discussed that this originates in the difference in the spreading width of the GDR 
between the two models. 
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FIGURES 
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FIG. 1. Strength functions of isovector dipole resonances in 18 calculated in TDDM (solid 
line) and in TDHF (dotted line). 



0.8 r- 
0.7 - 
0.6 - 




E [MeV] 



FIG. 2. Strength functions of isovector dipole resonances in 20 O calculated in TDDM (solid 
line) and in TDHF (dotted line). 
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FIG. 3. Strength functions of isovector dipole resonances in 22 calculated in TDDM (solid 
line) and in TDHF (dotted line). 
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FIG. 4. Strength functions of isovector dipole resonances in 24 calculated in TDDM (solid 
line) and in TDHF (dotted line). 
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FIG. 5. Isotope dependence of the low-lying El strength. Square, triangle, rhombus, and 
circle indicate the experiment 2 ), the shell-model calculation 3 ) , and the TDDM results with 
v = -230MeVfm 3 and -330MeVfm 3 , respectively. 
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